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Abstract 

\ Double-beta decay matrix elements (ME) for 76 Ge are calculated with dif- 

ferent quasi random phase approximation (QRPA)-based methods. First, the 
ME for the two-neutrino mode are computed using two choices for the single 
particle (s.p.) basis: i) 2 — 4huj full shells and ii) 3 — ATilu full shells. When 
calculated with the renormalized QRPA (RQRPA) and full-RQRPA their values 
are rather dependent on the size of the single particle basis used, while calcu- 
lated with proton-neutron QRPA (pnQRPA) and second-QRPA approaches such 
a dependence was found to be small. The Ikeda sum rule was well fulfilled within 
^ ■ pnQRPA for both choices of the s.p. basis and with a good approximation within 

second-QRPA, while the RQRPA and full-RQRPA methods give deviations up 
to 21%. Further, the ME for the neutrinoless mode are calculated with the pn- 
QRPA, RQRPA and full-RQRPA methods. They all give close results for the 
calculation with the smaller basis (i), while for the larger basis (ii), the results 
differ significantly either from one method to another or within the same method. 
Finally, using the most recent experimental limit for the Qv/3/3 decay half-life of 
76 Ge a critical discussion on the upper limits for the neutrino mass parameter 
obtained with different theoretical approaches is given. 

Pacs: 21.60.Jz (Hartree-Fock and random-phase approximations). 
Pacs: 23. 40. He (Relation with nuclear matrix elements and nuclear structure) 
Pacs: 23.40.Bw (Weak interaction and lepton (including neutrino) aspects) 
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1 Introduction 



Since the nuclei which undergo a (3(3 decay are generally rather far from the closed 
shells, the QRPA-based methods have been extensively employed for computing 
ME involved in the theoretical description of this process []i[|- p7| . Moreover, in 
spite of the recent progress of the shell- model and/or Monte-Carlo shell model 
techniques [^1] these methods also remain, at least for the next future, the only 
available for treating nuclear systems which are far away from the closed shells. 
The pnQRPA Q was the first adaptation of the standard QRPA for nuclear 
charge-changing processes. One of its most important achievements was after 
the pioneering work of || the success in explaining the suppression mechanism 
of the two-neutrino double beta (2v(3(3) decay ME reducing thus the large 

discrepancy existing until that moment between the theoretical and experimental 
(3(3 decay half-lives. However, this method faces the problem of a strong depen- 
dence of these ME on the renormalization of the particle-particle component of 
the residual interaction. Namely, if one represents the ME as function of the 
particle-particle interaction strength (usually denoted by g pp ), one observes that 
they decrease rapidly and cross through zero in a region of physical values of this 
constant, making the task of fixing it adequately difficult. To overcome this prob- 
lem several further developments of this method have been advanced during the 
recent past. We remind here: the appropiate treatment of the particle-number 
non-conservation ||10fl - |13|| , the inclusion of the proton- neutron pairing [O], the 
double commutator method [12], [14|, computation of the transitions to excited 
final states [14]-[|ll|, [24] as well as the development of approaches going beyond 
the quasi-boson approximation || , [|l6| , |l^] , |2(| , [E7| . At this point it is worth 
mentioning a nice feature of these higher-order QRPA approaches: the like- and 
unlike-nucleon residual interactions appear both in the next higher-order terms 
beyond pnQRPA, obtaining thus a more realistic picture of the competition be- 
tween them in producing a (3(3 decay. As a result, calculated with these methods 
the ME become more stable against g pp and the RPA break-down point is shifted 
towards the region of un-physical values of this constant. This is why, the fur- 
ther improvement of such approaches seems to be the most promissing line of 
development for treating the nuclear ME involved in the (3(3 decay process. 



The first method which has included higher-order terms beyond pnQRPA was 



developed in |J and further, applied with some modifications in [15|-|17[|, |22[ |. In 
this approach the extension of the pnQRPA was done using a boson expansion 
of both the phonon operators and transition /3 ± operators and retaining the next 
order in this expansion beyond the quasi-boson approximation (QBA). Also, this 
method allowed, for the first time, the computation of (5(5 decay rates to excited 
final states. An alternative approach for extending pnQRPA is based on the 
idea of replacing the uncorrelated QRPA ground state (g.s.) by a correlated 
g.s., in the calculation of the expectation value of the commutator of the two 
bifermion operators involved in the derivation of the QRPA equations. The 
expectation values of the number operator in the QRPA correlated g.s. are 
introduced in the quasi-boson commutators of the pair operators and this leads 
to a renormalization of the QRPA forward- and backward-going amplitudes. This 
method (called RQRPA) was first developed in refs. p3]-p3] for the standard 



QRPA and adapted later on for charge-changing processes in ref. [18]. Within 
the RQRPA a stabilization of the ME against g pp and a shift of the RPA break- 
down point towards larger (un-physical) values of this constant are also observed. 
However, this method has a main inconvenience consisting in an undesirable 
violation of the Ikeda sum rule (ISR) . Some refinements in the way of calculating 



the averages of the quasiparticle number operator are proposed [PQ|-p3l, but they 
result in a rather small reduction of the violation. 

In this paper we want to make a study of the (5(5 decay nuclear ME of 76 Ge 
calculated with different QRPA-based methods with the same set of parameters 
and for both two neutrino and neutrinoless modes. The motivation of such 
a study is given by some discrepancies concerning their values which are still 
found in the literature, where similar calculations have been performed. First, we 
calculated the nuclear ME involved in the 2v(5(5 decay mode using the pnQRPA, 
RQRPA, full-RQRPA and second-QRPA methods. One of our goals was to see to 
what extent the size of the single particle basis influences the values of these ME 
and how one can explain the differences between various calculations. A similar 



study has been made in [20] but only for the neutrinoless mode. Another point 
we have focused on was to check the Ikeda sum rule (ISR) in the framework of 
the above mentioned methods. Particularly, we would like to compare, under the 



same conditions of calculation (i.e. same parameters and s.p. basis), the various 
degrees of deviations obtained with these different methods and give possible 
explanations for the differences. Further, the ME for the neutrinoless mode are 
calculated with the pnQRPA, RQRPA and full-RQRPA methods. The results 
are found to be close to each other for all three methods in the case we used a 
smaller s.p. basis (9 levels), while for a larger one (12 levels) the results differ 
significantly either within the same method or from one method to another, for 
the two choices of the s.p. basis. Then, using the most recent experimental results 
for the two-neutrino and neutrinoless (3(3 decay half- lives of 76 Ge [p9| , we 
fixed first the g pp constant and then extracted new limits for the neutrino mass 
parameter. Finally, we give a critical discussion on the values of this parameter 
obtained with different theoretical methods. The paper is organized as follows: 
in section 2 we will give a short comparative description of the QRPA-based 
methods that we used for the calculation. The results are presented in section 3 
and the section 4 is devoted to the conclusions. 

2 Formalism 

In the QRPA-based methods one assumes the nuclear motion to be harmonic 
and the excitation QRPA operator may have the following general expression: 



r jM» = E {x™,(kJ,jnAl^(k^,J,M) + Y™,(kJ,r)A^(k,l,J,M) 

k,l,fjb<ti' 

(2.1) 

Here the summation is taken with k < I if /i = //. X m and Y m are the forward- 
and backward-going QRPA amplitudes and A, A^ the pair quasiparticle operators 
coupled to angular momentum J and projection M: 



l^(M,J,M)=Af(W//) £ C^ knmi a\ krn a\, lrni 



A^(k,l,J,M) = (-) J - M A^(k,l,J,-M) (2.2) 

N is a normalization constant, which is different from unity only in case when 
both quasiparticles are in the same shell [2C], fi,fj,' = 1,2 and 1 = protons, 2 = 



neutrons. Using the equation of motion method one can derive the pnQRPA 
equations which, in the matrix representation, may be written as: 



A B 
B A 



X r 



U 

-u 



x r 

Y" 



(2.3) 



j* \ / \ / J w 

where the matrices A, B and U have the following expressions: 



Aj(fik, vl;fi'k', v'l') = (0+ PA \ \A^(k,l,J,M),[H,Al, u ,(k',r,J,M)]\ \0+ pA ) 



Bj(^k, vl; fi k ; v'l') = (0+ | A, u (k, I, J, M), [Au> v >(k , I', J, M),H] 
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U = Kpa\ \A^(k,l,J,M),[Al v ,(k',l',J,M)] 
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RPAI 



RPAI 

(2.4) 
(2.5) 



Here the fiji are the QRPA excitation energies for the mode J n . 

Within the pnQRPA the QBA is assumed, i.e. the quasiparticle operators 
A, At are bosons and satisfy exactly the boson commutation relations: 



A^(k,l,J,M),Al, u ,(k',l',J,M) 
M(kfi, lv)M{k'^\ l'v') (s^S^SwSw - Wfe'fef) it+JrJ ) (2-6) 

In this way the Pauli principle is violated and this is a serious drawback of this 
method. To improve the situation in the RQRPA method the A, A' operators 



are renormalized [18], [20]: 



where the D^ kuk > matrices are defined as follows: 
D^kvk'jK = N{kn,lv)N(k! n' ,l'v') (5w'6 vv /8 kk/ Su> - 5^'5 v ^5i kl 5 kV (-) jk+jl 

By inspecting (2.6) and (2.8) one observes that by this renormalization one goes 
beyond the QBA by taking into account the next terms in the commutator 
relations of the A, A^ operators which are just, essentially, the proton and neutron 



(2.7) 



(2.8 



number operators. It is worth mentioning that they are taken into account within 
RQRPA only by their averages on the RPA g.s.. The renormalization of the A, 
operators is further carried onto the RPA amplitudes, on the A, B matrices and 
on the RPA phonon operator also obtaining a renormalization of them: 



X m = D l/2 x m . yrn = £,l/2ym . jjn = £-1/2^-1/2 . gm = Q-l/Zgp-l/l 

(2.9) 



T Jm«= E [xZ'(kJ,JnA^(kJ,J,M)+Y^(kJ,r)A^(k,l,J,M) 

(2.10) 

To calculate A and B we need to determine the renormalization matrices D. 
This is done by solving a system of non-linear equations for them by an iterative 
numerical procedure. As input values one can use their expressions in which the 
averages of the number operators are replaced by the back-forwarded amplitudes 
obtained as initial solutions of the QRPA equation. 

In QRPA- type methods, before starting the RPA procedure, we need the 
occupation amplitudes (u, v) and the quasiparticle energies, in order to get the 
image of the RPA operators in the quasiparticle representation. This is done by 
solving the HFB equations, which may include, in the general case, both like- 
and unlike-nucleon pairing. When one includes only like-nucleon pairing in these 



equations, the QRPA procedure described above was called RQRPA |jlq] , [20|, 



1 25], |p7| 1, while when both types of the pairing interaction are included it was 



called full-RQRPA f2§, J2§. On the other hand, if one takes the D = 1 we get 
back the QBA and these methods become pnQRPA and full-QRPA, respectively. 

In the second-QRPA method the principle of including higher-order correc- 
tions to the pnQRPA and restoring partially the Pauli principle is different. Here, 
the two quasiparticle and the quasiparticle-density dipole operators are expanded 
in a Beliaev-Zelevinski series 

4>«) = E K'°H(*0 + A^rt^k)) (2.ii) 



fclfc2 

(2.12) 



where 



EnJM t . 

jpm p j„m„ Ul i v m v 



^(pra) = (-) J ~ M B lfl (pn) (2.13) 

The boson expansion coefficients ^l^ 1 ' ), A^ , B^ , B^ are determined 
so that the equations (2.11)-(2.12) are also valid for the corresponding ME in the 
boson basis. 

Further, the transition (3^ operators in the quasiparticle representation can 
be expressed in terms of the dipole operators and 

13- (k) = e k A\^k) + d k M» + mB\^k) + f) k B^ 

#(*0 = - (M„(*0 + 9 k A llx + f] k Bl(k) + m B lfl ) (2.14) 

where 



6 k = 3jL(j p \\ a \\j n )U p V n - h = lL(j p \\a\\j n )U n V p - j = ^2j + l 



m = ■j=(J P \\<r\\jn)UpU n ; Vk = ^(j P \\cr\\jn)V p V N (2.15) 

Using the boson expansions (2.11)-(2.12), one also gets expressions of the transi- 
tion operators beyond the quasiboson approximation. Thus, in the second-QRPA 
method, the higher-order corrections to the pnQRPA are introduced not only in 
the RPA wave functions (by improving the phonon operator with additional cor- 
relations), but also in the expressions of the (3^ operators, and the procedure is 
now more consistent. The additional terms will have, of course, an influence on 
the ME calculation of these operators . 

Further, we give the factorized forms of the two-neutrino and neutrinoless (3(3 
decay half-lives that we used in our calculations : 



rp2v 

1 l/2 



F 2v \M% T \ 2 



(2.16) 



where F 2u is the lepton space phase and 

(0;||ar-||l+ fe )(l+|l+)(l+|| f 7r-||0+) 



E 

l,k 



Ei + Qpp/2 + m e -E 



(2.17) 



In (2.17) 1, k denote the two different sets of 1 + states in the odd-odd nucleus ob- 
tained with two separate RPA procedures applied onto the g.s. of the initial and 
final nuclei participating in the (5(3 decay. £/ is energy of the I — th intermediate 
1 + state, and Eq is the initial g.s. energy. 



rnOv 
1 l/2 



— c„ 



(2.18) 



where (m v ) is the effective neutrino mass and 



— p u 

mm — 1 1 



Mqt -[ — ) 



(M° 



(2.19) 



Ff" is the phase-space integral and Mqj, and Mp u are Gamow- Teller and Fermi 
matrix elements. 



3 Results 

3.1 Two-neutrino double beta decay 

First, we have performed calculations of the nuclear ME involved in the 2v(3(3 
decay mode of 76 Ge using the pnQRPA, RQRPA, full-RQRPA and second-QRPA 
methods. For the s.p. basis we used two choices. We included: i) the 12 
levels belonging to the full sd, pf and sdg shells, taking thus w O as core and 
ii) the 9 levels belonging to the full pf and sdg shells and taking thus 40 Ca as 
core. The single particle energies have been obtained by solving the Schrddinger 
equation with a Coulomb-corrected Woods-Saxon potential. For the residual 
two-body interaction there was taken the Brueckner G-matrix calculated from 
a Bonn-OBEP. The quasiparticle energies and the BCS occupation amplitudes 
were derived by solving the HFB equation without and with proton-neutron 
pairing, separately for the initial and final nuclei, with both choices of the s. p. 
basis. For a complete calculation we included in the model space the states with 
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1.138 


1.214 


1.653 
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1.031 


1.194 


1.502 



Table 1: Values of the renormalizing constants. 



all the multipolarities J n . The renormalization constants were chosen as follows: 
g pp = 1.0 for all the multipolarities, except the 1 + channel for which it was 
left as a free parameter, and g p h = 1.0 for all the multipolarities except the 2 + 
channel where it was fixed to 0.8, since for larger values the p-h interaction in this 
channel is too strong producing the collapse of the RPA procedure. The value of 
all the constants, including those which renormalize the pairing interactions are 
presented in the Table 1. 

In Fig. 1 we displayed the Mq^ (in MeV~ l ) as function of g pp calculated 
with pnQRPA and second-QRPA methods. The two curves for each method 
represent the calculations performed with the two different s.p. basis. In the 
figure is also drawn the line representing the ME value corresponding to the 
latest experimental 2u[3/3 decay half-life of the 76 Ge, obtained by the Heidelberg- 
Moscow experiment: Tf% = 1.55 x 10 21 yr (@]). 



As it was already observed in previous calculations [||, |TI| the point where QRPA 
breaks down is pushed to higher values of g pp in the framework of the second- 
QRPA method as compared with the pnQRPA. The calculation also shows that, 
within these two methods the values of the ME do not depend significantly on 
the size of the s.p. space, especially in the region around the experimental value. 
The values of g pp which fit the best this experimental value are: 0.94 for both 
calculations performed with pnQRPA and 0.99 and 1.01 for the calculation with 
12 and 9 levels, respectively performed with second-QRPA. 




In Fig. 2 are displayed the same ME but calculated with RQRPA and full- 
RQRPA methods, Contrary to the previous calculations, in this case the differ- 
ence between the results obtained with different choices of the s.p. space, within 
the same method, is rather large. Indeed at values of g pp where M"q T crosses the 
line representing the experimental result, the values for the ME, obtained with 
the same method, differ from each other by up to 40 — 50% when the two different 
choices of the s.p. basis are used. The values of g pp for the best fits with the 
experimental value for the ME are: 0.977; 0.982 in the case of the RQRPA and 
0.975; 1.012 in the case of the full-RQRPA for calculations including 9 and 12 
levels in the s.p. basis, respectively. The differences persist when the s.p. basis 
was enlarged to 21 states. This different behavior of the calculation, obtained 
with RQRPA and full-RQRPA on one side, and with pnQRPA and second-QRPA 
one the other side, in connection to the choices of the s.p. basis, reflects the sen- 
sitivity of the former methods in computing the M"q T ME. One possible source 
of this sensitivity might have its origin in the numerical computation. Indeed, 
the self-consistent iteration procedures for solving the full RQRPA equations, 
for all the multipolarities, are very time-consuming and rather slow converging 
and might affect the precision of the calculation. For more reliable calculations 
improved numerical techniques are in our opinion further required. 




Spp 



Figure 2: a = full - RQRPA(12) ; b = full - RQRPA(9) ; c = RQRPA(9) 
d = RQRPA{12) 



pnQRPA 


RQRPA 


full-RQRPA 


second-QRPA 


w Ge 0.23 0.26 


20.06 21.34 


17.68 17.29 


2.7 3.7 


w Se 0.41 0.48 


19.66 19.94 


17.12 16.91 


3.13 4.18 



Table 2: The numbers represent the deviations (in percents) from the ISR calculated 
within the specified methods. The first (second) numbers in the rows represent the 
calculation with a s.p. basis containing 12 (9) levels, respectively. 



On the other hand, there are some theoretical arguments which could explain 
the different results for the ME obtained with RQRPA-like methods as compared 
to the other two. We will discuss them later, after having discussed the ISR. 

The ISR 

5_ - S + = E m |<0+ ||/3-||l+)| 2 - S m |(0+ ||/3+||l+)| 2 (3.1) 

was checked out in the framework of the four methods. The results are presented 
in Table 2, where the percentages of deviation from the correct value for each 
method and choice of the basis are given. The first values in the row represent 
the calculations with a s.p. basis with 12 levels, while the second numbers refer 
to the same calculation, but with 9 levels. 

One can see, as expected, that within the pnQRPA the ISR is very well 



fulfilled, while within RQRPA and full-RQRPA the deviations are between 17 — 
21%. One can also observe that within second-QRPA the deviations from the 



ISR are rather small, confirming the result reported earlier in refs. [15], 16], but 
at that time calculated including only the 1 + channel and 9 levels in the s.p. 
basis. We should mention that in the present second-QRPA calculation we did 
not take into account the three boson states which may introduce undesirable 
spurious states in the QRPA space. Including such states one also gets deviations 
up to 17% from the ISR. Looking for some theoretical arguments for a possible 
explanation of the different extent to which the ISR is fulfilled within the RQRPA 
and second-QRPA methods, one finds that one reason could be the existence of 
some inconsistencies of the RQRPA related to the way of partial restoration of 
the Pauli principle. 
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Figure 3: a = RQRPA{9) ; b = pnQRPA(9) ; c = RQRPA(12) ; d = pnQRPA(12) 

Indeed, as we already mentioned in section 2, within the RQRPA method 
the Pauli principle is partially restored for the operators A, A\ by taking into 
account the averages of the quasiparticle-number operators in their commutator 
relations. However, there is no justification to neglect them in the B, oper- 



ator commutation relations. Within the second-QRPA higher order corrections 
beyond pnQRPA are taken into account in the expression of these operators and 
moreover, such corrections are also introduced in the expressions of the op- 
erators. The effect of the additional terms combined with a larger boson space 
(in the second-QRPA the boson space enlarges from one to two boson states) 
reflects in a positive contribution to the ISR. On the other hand it is known that 
RQRPA underestimates the ISR, so this could be one possible explanation why 
the ISR is better fulfilled within second-QRPA. Another possible shortcoming 
of the RQRPA method is a lack of consistency between BCS and QRPA levels. 
While in the BCS still one assumes the g.s. to be the quasiparticle vacuum at 
the level of RQRPA we are dealing with the non-vanishing quasiparticle con- 
tent of the g.s. due to the additional scattering terms taken into account in the 
commutation relations [27|. 



3.2 Neutrinoless double-beta decay 

Further, we have performed a calculation of the neutrinoless ME using pnQRPA, 
RQRPA and full-RQRPA methods, also for the two s.p. basis. The M 0u as 
function of g pp calculated with the pnQRPA and RQRPA are displayed in Fig. 

3, while the same ME but calculated with the full-RQRPA are displayed in Fig. 

4. One observes that all the three methods used for calculation give different 
values of the ME for different choices of the s.p. basis. The differences between 
ME values calculated with 9 and 12 levels included in the s.p. basis, within the 
pnQRPA and RQRPA methods, are given by factors of 3 and 2.5, respectively, 
while for the full-RQRPA method the difference between the two calculations 
reduces to a factor of about 1.6. One also observes, that the values of the 
ME obtained with the three methods are close to each other (3.9-4.1) in the 
calculation with the smaller basis. When enlarging the basis to 21 levels, the 
result is close to that obtained with the basis with 12 levels. This again reveals 
the sensitivity of the RQRPA-type methods to the choice of the s.p. basis and 
seems to indicate a possible stabilization of the results for larger basis. However, 
a general conclusion about which basis is better to choose is difficult to give 
until we have not the whole image of a similar study performed on several other 
double-beta emitters. Another still open question is what are the results when a 



similar calculation is performed with the second-QRPA. 
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Figure 4: a = full - RQRPA(9) ; b = full - RQRPA(12) 



Finally, using the value of the g pp constant, fixed from the M"q T calculation 
to fit the most recent half-life, i.e. 1.0 (very close to the average value between 
the two second-QRPA and full-RQRPA calculations), and using the most recent 
experimental limit of the neutrinoless mode half-life for the Ge case (i.e. > 
1.9 x 10 25 yr (90% C.L.) reported by the Heidelberg-Moscow experiment pOfl ) 
we extract new upper limits for the neutrino mass parameter within the full- 
RQRPA method. We obtain the values = 0.407 eV and 0.625 eV, if we use in the 
calculation a s.p. basis with 9 and 12 levels, respectively. In addition, in table 3, 
besides our values for the ME and neutrino mass parameter, we present results 
of other calculations found in the literature. For a direct comparison between 
various results we use non-dimenssional values for all the ME taken from the 
references indicated in the table. Further, using the same phase space factor 

wc 



= 6.31 x 10 15 yr 1 |26|| and the same half-life reported in Ref. 



extracted upper limits for < m„ > corresponding to all values of the ME, in order 
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present work 




4.25 


4.26 


4.85 


1.57 


1.92 


2.80 


2.36(12) 3.62(9) 


< m u > [eV] 90% 


0.345 


0.328 


0.304 


0.940 


0.768 


0.527 


0.625 0.407 


<m u > [eV] 68% 


0.27 


0.26 


0.24 


0.73 


0.60 


0.41 


0.490 0.320 



Table 3: Neutrinoless M.E. and upper limits for the neutrino mass parameter for 76 Ge, 
calculated with a phase space F?" = 6.31 x 1CT 15 y~ l (Q2§|) and Ty 2 > 1.9 x 10 25 y 
(C.L. 90%), and 3.1 ■ 10 25 y (C.L. 68%) |IJ. The non-dimensional values of various 
ME are taken and reconverted, when necessary, from the indicated references. For the 
present work two values, representing the calculation with 12 and 9 levels for the s.p. 
basis, are displayed. 

to have a more complete image of the evolution of the theoretical predictions. 
Performed with pnQRPA, the older results are very similar, although they were 
calculated by different groups and with different numerical codes and parameters 
( (8| , (ll| ) . One also observes that their values are larger by a factor of about two 
than the values obtained by using the recent extensions of pnQRPA, RQRPA and 
full-RQRPA. However, it should be kept in mind that these last approaches do 
not fulfill the ISR and thus lead apriori to too small M.E i.e. too large neutrino 
mass limits . Our calculations performed with the larger s.p. basis give M.E. 
rather close to those of Ref. [f3^] were corrections due to the nucleon currents 
such as weak magnetism and pseudoscalar couplings to the amplitude of Qv(5(5 
have been taken into account. The use of the smaller s.p. basis yields a value for 



the M.E. which close to the earlier approaches ||, ||, [11|. It should be pointed 
out that corrections due to nucleonic currents mentioned above were not make 
in the present study. 

On the other side there are the values of the M.E. calculated with the shell- 
model in Refs. 0, [28| which differ from each other by a factor of 3. However, 



in our opinion, these calculations performed with the shell model are not, at 
present, reliable enough. This has been stressed also by pq ]. In ref. [Q] the 
calculations were performed with a rather crude shell-model code in a weak 
coupling approximation, at the computer performances of that time. Also, in 



calculations of ref. [28] some important orbits are missing, like some spin orbit 
partners, resulting in a violation of the Ikeda sum rule of about 50%, which 
means it should be expected that they give too small M.E.. 



4 Conclusions 



We have performed a calculation of the two- and zero-neutrino (3(3 decay matrix 
elements for the case of 76 Ge with the pnQRPA, second-QRPA, RQRPA and 
full-RQRPA methods, using two different choices of the s.p. basis. We can 
summarize the main results as follows: 

i) for the M"q T we got a significant dependence of the results on the size of the 
s.p. basis, in the case of the RQRPA and full-RQRPA methods, while the results 
obtained with pnQRPA and second-QRPA do not display such a dependence. 

ii) for the neutrinoless decay mode all the three methods used for the calcu- 
lation, i.e. pnQRPA, RQRPA and full-RQRPA, give differences between 9 and 
12 level calculations by factors 1.6-3. The values of the ME obtained with the 
three methods are close to each other for the calculation with the smaller basis, 
while they differ significantly when the calculation is done with the larger basis, 
i) and ii) reveal a sensitivity of the RQRPA methods to the size of the s.p. basis 
which is used. This could have its root in the numerical double-iteration proce- 
dure used in RQRPA-type calculations and in our opinion further improvements 
should be done in this respect. 

iii) we also check the ISR within the four methods and found it to be fulfilled 
with a good approximation within second-QRPA method, while with RQRPA 
and full-RQRPA the deviations are up to 21%. We found that this result is 
not much dependent on the size of the s.p. basis used. This result, besides the 
numerical arguments mentioned above, might also be explained by theoretical 
arguments related to the way the partial restoration of the Pauli principle is 
done within RQRPA. The restoration is made in the commutator relations of the 
operators A, by taking into account the averages of the quasiparticle-number 
operators in their commutator relations. However, there is no justification to 
neglect them in the B, operator commutation relations. However, this is done 
within the second-QRPA and moreover, in this method the next higher-order 
corrections beyond pnQRPA are also taken into account for the (3 ± operators. 
The additional terms give a positive contribution to the ISR, while as it is known 
RQRPA underestimates the ISR. 

iv) using the most recent reported neutrinoless half- life limit, and using the 



value of g pp fixed for the 2v neutrino mode calculation, we extracted the following 
new upper limits for the neutrino mass parameter. A critical comparison between 
various values of the M.E. found in the literature was performed. One observes 
a tendency of reducing these values according to the most recent calculations 
performed with RQRPA and full-RQRPA. One may conclude that the values of 
the M.E. involved in Ou/3/3 decay of 76 Ge can be reliably predicted within a factor 
of two. 

Finally we would like to stress that considering the various approximations 
made in the different calculations (violation of Ikeda sum rule in |2^] , pl[| , |35| , 
and this work (by 50, 20%)), and neglection of weak magnetism and pseudoscalar 
coupling in all approaches except in [|35| (another 30%), the tendency goes to a 
variation within the different approaches of only a factor of 1.5, and to clearly 
favouring the smaller deduced neutrino mass values. The < m v > values expected 
from 76 Ge decay would lie around 0.2 eV (68% C.L.) after the corresponding 
estimated corrections. 
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